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ABSTRACT 

We present theoretically-established values of the CO-to-H 2 and C-to-H 2 conversion 
factors that may be used to estimate the gas masses of external galaxies. We consider 
four distinct galaxy types, represented by M51, NGC6946, M82 and SMC N27. The 
physical parameters that best represent the conditions within the molecular clouds 
in each of the galaxy types are estimated using a \ 2 analysis of several observed 
atomic fine structure and CO rotational lines. This analysis is explored over a wide 
range of density, radiation field, extinction, and other relevant parameters. Using these 
estimated physical conditions in methods that we have previously established, CO-to- 
H2 conversion factors are then computed for CO transitions up to J = 9 — > 8. For the 
conventional CO(1-0) transition, the computed conversion factor varies significantly 
below and above the canonical value for the Milky Way in the four galaxy types 
considered. Since atomic carbon emission is now frequently used as a probe of external 
galaxies, we also present, for the first time, the C-to-Eb conversion factor for this 
emission in the four galaxy types considered. 
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1 INTRODUCTION 

The CO-to-H 2 conversion factor, X C o = AT(H 2 )/W(CO) 



cm (K km s 



where N(J&2) is the column density of 



molecular hydrogen (in cm -2 ) and TU(CO) is the velocity- 
integrated intensity of the CO J = 1 — ► rotational 
line (in Kkms -1 ), is widely used to estimate cloud masses 
within the Galaxy and in external galaxies. Howev er, as 
demonstrated in a previous study (|Bell et all 120061 . here- 
after Paper I), the value of Xco is highly sensitive to the 
local physical conditions, and may vary considerably even 
within the Galaxy. The appropriate value to use in ex- 
ternal galaxies has not so far received theoretical atten- 
tion. Observationally, studies of external galaxies have in- 
ferred values of Xco using a variety of techniques, including 
virial mass estimates when individual clouds can be resolved 
l|Brouillet et al.lll998h , dust continuum observations coupled 
wit h some assumptions regarding the gas-to-dust mass ra- 
tio (|Boselli. Lequeux fc Gavazzil 12003) and the use of neu- 
tral carbon emi ssion as a tracer when its a bundance can 
be determined (jlsrael. Tilanus fc Baai[2006). Such studies 
confirm the variable nature o f Xgo, which is seen t o dif- 
fer from galaxy to galaxy (e.g. iMalonev fc "B lack 198c]). The 
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ability to resolve molecular emission from individual clouds 
within nearby galaxies through interferometric techniques 
has created new opportunities for studyi ng the interstel- 
lar medium in different environments (e.g. ISakamoto et al.l 
1999-, iHelfer et al.ll20o3 ): however, many properties still de- 
pend critically on the masses derived for these clouds. The 
CO-to-Eh conversion factor therefore has an increasingly im- 
portant role to play in advancing our understanding of the 
properties of interstellar gas in external galaxies. Atomic 
carbon emission is now also routinely observed in external 
galaxies and has been pr oposed as an alternative tracer o f 
molecular mass (see, e.g.. iPapadopoulos. Thi fc Vitill2004f ). 
with a C-to-H2 conversion factor of the form X01 = 



AT(H 2 )/W(Ci) 



(Kkmi 



where W(Cl) is the 



velocity-integrated intensity of the [C 1] 609 \im fine struc- 
ture line (in Kkrns" 1 ). 

The purpose of this paper is to determine theoreti cally 
from fundamental arguments (as described in iPaper j) the 
appropriate values of Xco and Xci to use for each of sev- 
eral distinct galaxy types in which the physical conditions 
are significantly different from each other. We select a sam- 
ple of nearby galaxies that are representative of the various 
commonly observed morphological types and determine the 
physical properties that best describe the neutral interstel- 
lar gas within them by fitting observed line intensity ratios 
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for a range of far-infrared (FIR) and sub millimetre diagnos- 
tic lines. Using the methods of |Paper"H we then calculate 
Xco for various rotational transitions of the CO molecule 
and Xc i for the 609 fi.m fine structure transition of atomic 
carbon for each of the galaxies selected, taking the physical 
parameters from the best fitting model in each case. 

In Section [2] we discuss the sample of galaxies we have 
selected to study and introduce them individually. Section [3] 
describes our fitting procedure to determine appropriate 
physical parameters for the interstellar medium (ISM) in 
each of our sample galaxies. We present the best fit model 
for each galaxy in Section [3] and discuss the derived physi- 
cal parameters and the quality of the fits. In Section [5] we 
use these models to compute suitable CO-to-Hb and C-to- 
H2 conversion factors for each galaxy type considered and 
we summarize our conclusions in Section [6] 



2 THE GALAXY SAMPLE 

To produce CO-to-Eb and C-to-Eb conversion factors that 
are applicable to as wide a variety of galaxies as possi- 
ble, we select four nearby examples that are considered to 
be prototypical of the main galaxy types, possessing inter- 
stellar clouds that span a range of physical conditions and 
so are likely to require conversion factors that are some- 
what different to the canonical va lue for the Milky Way , 
Xco = 2x 10 20 cm" 2 (Kkms-V 1 fetrong fc Mattoxll 19961 ; 
iDame. Hartmann fc Thad dcus 200l]). The examples we con- 
sider are the grand-design spiral galaxy M51, the late-type 
spiral galaxy NGC 6946, which has a nuclear starburst, the 
irregular starburst galaxy M82 and nebula N27 within the 
Small Magellanic Cloud, a dwarf irregular galaxy. Table [1] 
lists the properties of our ch osen galaxies; their morphol ogies 
(from the RC3 catalogue of lde Vaucouleurs et al.|[l99ll ). the 
coordinates of their central positions (J2000 epoch), their 
radial velocities relative to the local standard of rest (LSR) , 
estimated distances (in Mpc) and the corresponding angular 
scale at that distance (in pc/arcsec). 

We restrict our study to the central positions of these 
galaxies (with the exception of the Small Magellanic Cloud, 
see below). This is partly due to the limited amount of rel- 
evant observational data in the required transitions. Most 
studies of the neutral ISM within external galaxies focus 
on single-beam observations or mapping of the brightest re- 
gions, typically their centres, which are known to exhibit 
strong molecular concentrations. These studies are further 
restricted due to the large telescope beam sizes involved 
(typically >10 arcsec, although interferometric mapping can 
reach angular resolutions of <~1 arcsec); these beams can en- 
compass whole ensembles of clouds that exhibit a range of 
properties. Determination of the physical conditions within 
these regions is then necessarily limited to deriving aver- 
age properties for an ensemble of clouds. In this section, the 
four galaxies in our sample are introduced and observational 
studies of the interstellar gas within them are summarized. 



2.1 M51 

The grand-design spiral of M51 (NGC 5194) appears nearly 
face-on (i ~ 20°) at a distance of 9.7 Mpc (see Table [1]), 
and is seen to be interacting with its smaller companion, 



NGC 5195, which lies 4.5 arcmin to the north. It has a metal- 
licity slightly above Solar, d etermined near its centre to be 
12+ l og (O/H) = 9.28±0.11 (|Zaritskv, Kennicutt fc Huchral 
1994). M51 has been c l assed a Seyfert 2 galaxy by 



Ho. Filippenko fc SargentI (|l997l '>. with a central AGN sur- 
rounded by a ~100 pc disc of warm and dense gas 
l|Matsushita et al] 1 19981 ). This galaxy has been studied ex- 
tensively and a large amount of observational data are avail- 
able, including an extended map of [Cn] 158 urn emission, 
taken with the Kuiper Airborne Observatory (INikola et al.l 
2001) and additional observations of FIR fine-structure lines 
toward the galaxy centre (|Negishi et al.l [2001) using the 
Long Wavelength Spectrometer (LWS) onboard the Infrared 
Space Observatory (ISO). Ground-based single-dish stud- 
ies mapping the low- lying transitions of 12 CO and 13 CO 
sho w the CO to be tightly confined to the spiral arms 
fe.g.lGarcfa-Burillo. Guelin fc Cernicharolll99 3; Nak ai et al.l 
1199411 . Maps of higher CO transitions. J = 3^2 and 4 -» 
3, have be en obtained with th e Hei nrich Hertz T el escop e 
(HHT) by iNieten et al] (|l999l ) and iDumke et ail i|200ll ). 
showing warm molecular gas to be present out to a galac- 
tocentric radius of 5 kpc or more within these arms. Re- 
cently, the [C 1] 609 um fine structure line has also been ob- 
served towards t he galaxy centre and several positions along 
the spiral arms (iGerin fc Phillipsll200o1 ; llsrael fc Baasll2002l ; 
iKramer et al.ll2005l ). which show that the neutral carbon 
emission is also extended. Apertu re syn thesis maps w e re ob - 
tained by ISakamoto et alj (| 19991) and iRegan et ail (|200ll ) 
at resol utions of 4-6 arcsec in t he CO (1-0) line and, more 
recently. [Matsushit a et aU (|2004T ) have mapped the inner re- 
gion in CO(3-2). 



2.2 NGC 6946 

This nearby low-inclination (i ~ 30°) galaxy, classified as 
a spiral with a weak bar (SAB) in the RC3 catalogue 
(jde Vaucouleurs et al.lll99lh . is at a distance of 5.5 Mpc (see 
Table [TJ) and has an abund ance gradient rilog ( Q/H)/ri.R = 
-0.089 ± 0.003 dexkpc -1 |Bellev fc Rovlll992) . similar to 



that of the Milky Way, and a metallicity of 12 + lo g (O/H) 
9.13 ± 0.15 near its centre ilZaritskv et al.lll994l). I t has 



nuclear starburst (see, e.g.. lEngelbracht et al.lll996l ) and a 



small inner bar, seen clearl y in recent interferomet ric maps 
of CO(1-0) and CO(2-l) (|Schinnerer et al.ll2006l ). that is 
driving an inflow of molecular gas into the central region 
of intense star formation. The galaxy shows a pronounced 
open spiral pattern, the north-eastern arm being the bright- 
est at all wavelengths, indicating a site of vigorous ongoing 
star formatio n. One of the first galaxies to be mapped in 
CO emission (|Morris fc Lolll978T l. NGC 6946 has been rel- 
atively well studied in the lo wer CO transitions, both to- 
ward s its central region (e.g. Weliachcw, Casoli & C ombesl 
1988) and towards its spiral arms 1 Casoli et al.lll99"oT ). More 
recently, maps of the central region have been made in the 
J = 2 — > 1, 3— > 2 and 4 — > 3 transitions of CO and in the 
[C i] 609 urn line l|Mauersberger et al.1 fl999l ; llsrael fc Baasl 
2001). All these observations show that the gas is predom- 
inantly molecular within 5 kpc of the centre and that CO 
emission is rather extended over the galaxy, but is most con- 
centrated in the central region. 
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Table 1. Properties of the four galaxies adopted for this study. 



Galaxy 


Morphology" 


R.A.(J2000) 


Dec.(J2000) 


VhSR 


Distance 


Angular Scale 






(hh mm ss) 


(° ' ") 


(kms" 1 ) 


(Mpc) 


(pc/arcsec) 


M51 


SA(s)bc pec 


13 29 52.7 


+47 11 42.0 


464 


9.7" 


46.5 


NGC 6946 


SAB(rs)cd 


20 34 51.2 


+60 09 17.5 


55 


5.5 C 


26.7 


M82 


10; Sbrst 


09 55 52.3 


+69 40 45.9 


175 


3.2 d 


15.5 


SMC N27 


dSB(s)m pec 


00 48 21.1 


-73 05 29.0 


115 


63 kpc e 


0.3 



References: a lde Vaucouleurs et al 
d lDumke et alj l l200ll ); e lRubio et al 



" ISandage fc Tammannl dl975h : c lTultvl l|l988l l: 



2.3 M82 

The irregular starburst galaxy M82 (NGC 3034) is per- 
haps the most studied of all nearby galaxies due to its 
spectacular display of active star formation, which domi- 
nates its evolution and is thought to have been triggered 
about 100 Myr ago by tidal interaction with its companion 
galaxy M81. Originating in the nucleus, the starburst ac- 
tivity within M82 is currently seen to be propagating into 
the molecular rings, disrupting the surrounding ISM and 
causing the insterstellar clouds to fragment. The gravita- 
tional interaction is also likely to hav e produced the ste l- 
lar bar seen in the central ~400 pc |Telesco et al.|[T99ll ). 
As a result of its vigorous star formation, M82 is also one 
of the most luminous nearby galaxies in the infrared, with 
Lxr ~ 4 x 10 10 I/©. The galaxy has been the subject of 
many observational studies, including extensive mapping 
in both low and high (up to J = 7 — » 6) transitions 
of GO (e.g. ISutton. Phillips fc Massonl Il983l; IWhite et al " 



1994 1; iMao et alj l2000l ; ISeaauist. Lee fc Moriartv-Schieven 
2006), comprehensive FIR spectroscopy IjColbert et al 



1999) and mapping in a v ariety of molecular speci es, in- 
cluding HCO+ and HCN dSeaauist fc Fraverl |200Ch, SiO 
|Garcfa-Burillo et alJlioOlh and CO+ (|Fuente et al.ll2006l ). 
These observations demonstrate the presence of PDR chem- 
istry over large regions of the galaxy centre and strong, ex- 
tended CO emission arising from a warm and diffuse ISM. 



2.4 Nebula N27 in the Small Magellanic Cloud 

At just under 63 kpc away, the Small Magellanic Cloud 
(SMC; NGC 292) is one of the nearest examples of dwarf ir- 
regular galaxy, presenting an excellent opportunity to study 
the properties of interstellar gas under conditions very dif- 
ferent to those in the Milky Way. Its metal content is 
known to be particularly low, wi th an oxygen abundance 
of 12 + log [O/H] = 7.96 ± 0.11 l| Vermeil fc van der Hulst] 
2002) implying a metallicity of about one tenth Solar, and 
it possesses a gas-rich interstellar medium that is subjected 
to strong FUV fluxes due to t he presence of a la rge popu- 
lation of young massive stars l|Westerlundl Il99p| ) . As part 
of a three-galaxy system that includes the Large Magel- 
lanic Cloud (LMC) and the Milky Way Galaxy, the SMC 
has experienced a disruptive history; an encounter with the 
LMC some 200 Myr ago is likely to b e responsible for the 
curre nt burst of star-forming activity dGardiner fc Noguchil 
1996). Such rapid star formation in a low metallicity envi- 
ronment is also thought to be representative of the condi- 
tions within primordial galaxies. In the absence of an obvious 
central feature, we have chosen to adopt the well-studied 
region surrounding the N27 nebula as our example of the 



typical environment within the SMC. The N27 nebula is 
a bright and relatively compact Hll region located in the 
south-western bar of the SM C and is associated with the 
IRAS point source LIRS 49 (|Schwering fc Israeli [19891 ). It 
has bee n mapped in the J = 1 — > a n d J — 2^1 lines of 
CO by iRubio. Lequeux fc Boulangerl (|l993f ) and observed 
in the J = 3 — > 2 CO line and transitions of CS, HCO + 
and HCN bv lHeikkila. Johansson fc Olofssor] (ll999T ). More 
recently, the [C 1] 609 urn line has been d etected in this re- 
gion for the first time (|Bolatto et al.ll2~000l ) and the FIR fine 
structure lines have been observed using the ISO LWS. 



3 PHYSICAL PARAMETERS FOR THE 
GALAXIES 

The atomic fine structure and CO rotational emission lines 
serve as excellent diagnostics of the physical conditions 
within PDRs, since they are the dominant cooling lines 
within the warm gas at the surfaces of molecular clouds. 
Their range of critical densities and energy level spacings 
mean that, when used in combination, the intensities of these 
lines can be used to constrain the density, in cident FUV flux 
and s urface temperature of the gas (see, e.g.. lKaufman et al.l 
I1999T) . 

In order to derive appropriate physical parameters for 
the galaxies in our sample, we compare emission line inten- 
sity ratios predicted by PDR models to those detected in 
each of the four example galaxies. The best fitting models 
are determined by performing a \ 2 goodness of fit analy- 
sis of the predicted and observed line ratios for a grid of 
PDR models covering a large region of parameter space. 
Section T3 . 1 1 describes the grid of models used for comparison 
with the galaxies in our sample, the calculation of appro- 
priate line ratios from the observed line intensities and the 
calculation of \ 2 values, from which the best fit models are 
determined. The observational data obtained from the liter- 
ature and the resulting line ratios adopted for each galaxy 
are discussed and the corresponding best fit models are pre- 
sented in Section [4] These models are then used in Section [5] 
to derive suitable CO-to-H"2 and C-to-Eb conversion factors, 
Xco and Xci, for each galaxy type considered. 



3.1 Best fit PDR models using x 2 analysis 

The short computation times required by the UCL_pdr code 
make it ideally suited to the calculation of large grids of 
models covering substantial regions of parameter space. The 
models predict the integrated intensities of the fine structure 
and CO rotational lines emitted by a PDR with the adopted 
physical parameters. Each model produces line intensities 
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Table 2. The region of parameter space covered by the 
grid of PDR models. 



Parameter Range of values 



Cloud density (cm -3 ) 


10 2 < 


n ^10 b 


Incident FUV flux (Habing) 


10 1 s: 


Go < io 5 


Age of the cloud (yr) 


io 4 < 


t ^10 8 


Cloud size (mag) 


OsC 


A v sCIO 


Metallicity 


Z/Zq 


= 1.0,0.5,0.1 


C.R. ionization rate (s — 1 ) 


C = 5x,50x,500xl0- 17 


Turbulent velocity (kms -1 ) 


Vtuxb z 


= 1.5 



as a function of total cloud depth (at visual extinctions of 
< A v < 10 mag) for cloud ages of 10 4 , 10 5 , 10 6 , 10 7 and 
10 8 yr. The standard version of the ucl_pdr code was used 
to calculate the grid of models to which the y 2 analy sis is 
then applied. The code is described in detail in lPaper 1 The 
gas is initially assumed to be fully atomic in the models, with 
all carbon in C + , and the evolving chemistry and thermal 
balance are then computed for subsequent time-steps. The 
PDR is treated as a one-dimensional semi-infinite slab of 
uniform density. We adopt the standard chemical network 
used in the ucl_pdr code, neglecting freeze-out of atoms 
and molecules on to grains. The metallicity-dependence of 
the gas-phase elemental abundances, dust grain and PAH 
number densities, photoelectric heatin g rate an d H2 forma- 
tion rate on grains are as described in lPaper II Table [2] lists 
the range of parameter values covered by the grid of models. 
Homogeneous clouds with hydrogen nuclei number densities 
of 10 2 < n < 10 5 cm" 3 , subject to incident FUV fluxes 
of 10 sC Go < 10 s times the standard lHabind (|l968l ) ISRF 
have been computed at metallicities Z/Zq of 1, 0.5 and 0.1, 
and cosmic-ray ionization rates £ of 5xl0" 17 , 5xl0" 16 and 
5xl0" 15 s" 1 . A turbulent velocity of i> tU rb = 1.5 kms" 1 
is adopted in all the models, typical of individual Galactic 
molecular clouds (e.g. iHollenbach fc TielensHl999l ). We do 
not consider the variation of this parameter in our attempts 
to fit the observed emission from each galaxy, since the line 
ratios are believed to b e fairly insensitive to small changes i n 
the turbulent velocity l|Wolfire. Hollenbach fc Tielenslll989h . 
Furthermore, few estimates exist for the degree of turbulent 
motion within molecular clouds of other galaxies, since ob- 
servations of extragalactic emission lines do not generally 
resolve the velocities of individual clouds. 

Comparing PDR model predictions to extragalactic ob- 
servations is never straightforward. Individual molecular 
clouds are rarely resolved in single-dish observations, partic- 
ularly so in the case of FIR fine structure lines, which have 
typically been observed with the ISO LWS with a beam size 
of ~ 70-80 arcsec. Several phases of the ISM are therefore 
contained within a single beam and the line intensity ra- 
tios will represent the combined emission from these phases. 
The physical properties derived from comparison with PDR 
models are then the average properties for the molecular 
clouds within the region, although they may be biased to- 
wards those of the brighter components. By using line ratios, 
rather than fitting the line intensities directly, the beam area 
filling factors of the two emission lines are cancelled out, as- 
suming that they come from the same clouds. 

In our fitting of the observed fine structure and CO 
rotational lines, we will consider the following line intensity 
ratios: 



(1) [Ol] 63um/[Cn] PDR 

(2) [Or] 145 urn/ [Or] 63 \im 

(3) [Cii] PD r/[Ci] 609 urn 

(4) [Cii] PD r/CO(3-2) 

(5) [C 1] 609 Ltm/CO(3-2) 

(6) CO(3-2)/CO(l-0) 

when the necessary line intensities are available in the liter- 
ature. The calculation of [C ii]pdr is described below. 

In order to derive meaningful ratios, the various line 
intensities must share a common spatial resolution. Given 
that the lines are observed with different beam sizes, it is 
necessary to convolve the integrated intensities to the spatial 
resolution appropriate to the largest beam size, typically 
that of the ISO LWS, which has an average half-power beam 
width (HPBW) of ~75 arcsec. When observational data at 
the necessary resolution are not available in the literature, 
we perform the convolution using the scaling factor: 

where #t>, > and 9h,f are the initial (observed) and final beam 
sizes (HPBW in arcsec), and 6 B , X and StV are the source 
sizes along the major and minor axes (FWHM in arcsec). 
Of course, if the source is circular then these two axes will 
have the same dimensions. The observed intensity is then 
multiplied by this factor to obtain the intensity at the de- 
sired resolution. 

Since [C 11] 158 u.m emission arises from both ionized 
and neutral gas, the fraction originating from the ionized 
medium must be subtracted before a comparison with PDR 
models can be made. We have chosen to estimate the frac- 
tion of [C 11] emission from the ionized medium using the 
observed [Nil] 1 22 urn fine structure line, in the manner of 
iMalhotra et~atl l|200ll ) and lContursi et all (|2002m . The [N 11] 
emission occurs exclusively in the ionized gas, so by assum- 
ing a constant ratio between the [C 11] and [N 11] line in- 
tensities in the diffuse ionized medium (DIM), taken to be 
the dominant ionized component, the contribution of ionized 
gas to [C 11] emission can be estimated. The [C n]/[N 11] ratio 
derived by these authors is 4.3, based upon observations of 
ionized gas in the Milky Way and using Solar abundances for 
carbon and nitrogen. In reality, this ratio is likely to vary 
from region to region and the metallicities of the galaxies 
studied here are known to differ. However, the abundance ra- 
tio has been shown to be indep endent of metallicity in both 
normal and irregular galaxies l|Garnett et all 1 19991) and in 
view of the lack of more reliable estimates of the ionized gas 
contribution, the use of this ratio is justified. The intensity 
of [C 11] emission originating from PDRs is then 

7 [Cn] PDR = I[Cn] ~ 4.3 x J [Nn] [ergs -1 cm" 2 sr" 1 ], (2) 

where 7[c „j and 7[ N n i are the observed intensities of the [C 11] 
158 [im and [N 11] 122 [im lines. 

Finally, the [C 1] and CO line intensities are usually 
measured as velocity-integrated antenna temperatures (in 
Kkms -1 ). These values must therefore be converted to in- 
tegrated intensities (ergs -1 cm -2 sr -1 ) before they can be 
compared to the FIR fine structure lines and model predic- 
tions. This is achieved using the relation: 

I=—j- I Ta<1v [ergs -1 cm -2 sr" 1 ], (3) 
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where v is the transition frequency (in Hz) and f T\ dv is the 
antenna temperature integrated over velocity (in Kkms -1 ). 

Once the relevant line intensity ratios have been deter- 
mined for a particular galaxy, we calculate the reduced \ 2 
goodness of fit for each model in the grid, 



2 

X = 



{No 



rjobs r> m 
Hi — It, 



(4) 



where A^bs is the number of independent line ratios com- 
pared to observations, iV par is the number of free parameters 
in the models and, for each line ratio i considered, R° hB is 
the observed value, Rf lod is the model value and m is the 
estimated error on the observed value. We take the typical 
observational error to be ~20 per cent of the line inten- 
sity, including calibration uncertainties, leading to a com- 
bined error of ~30 per cent on the line ratios. In addition 
to these uncertainties, errors introduced by convolving the 
intensities to the required resolution and by correcting for 
the ionized gas contribution also play a significant role in de- 
termining the value of \ 2 ■ These errors are hard to quantify 
and have therefore been neglected here, but their inclusion 
would serve to generally decrease the value of \ 2 obtained 
for each model. Uncertainties in the contribution of the ion- 
ized gas emission would only affect line ratios involving the 
[C n] 158 vim line, reducing the significance of these ratios 
in constraining the best fit model. Note that three of the 
line ratios are linked and so not independent, [C ii]pdr/[C i] 
609 nm, [Ci] 609 u.m/CO(3-2) and [C ii] PD r/CO(3-2), thus 
reducing the number of independent line ratios JV b s by one. 

The best fit model for the physical conditions in each 
galaxy is therefore the model with the minimum reduced \ 2 
value. In the next section we list the observed and predicted 
line ratios for each galaxy in our sample and present contour 
maps of reduced \ 2 across the grid of models for each galaxy. 
The physical parameters derived from the best fit models are 
discussed for each of the galaxies in our sample. 



4 ADOPTED LINE RATIOS AND BEST FIT 
MODEL PARAMETERS 

The line intensities used to derive the relevant line ratios 
are taken from the literature, applying the convolution for- 
mula in equation [1] when necessary in order to obtain a set 
of integrated intensities at a common resolution for each 
source. All the FIR data included have been acquired using 
the I SO LWS, with a typical HPBW of 70-80 arcsec i|Llovdl 
2003), corresponding to spatial scales of ~l-3.5 kpc at the 
distances of these nearby galaxies. All other observations at 
higher spatial resolution have therefore been convolved to a 
beam size of 75 arcsec (fib = 1-1 x 10~ 7 sr). 

Table [3] lists the adopted line intensity ratios for each 
galaxy, to which the predicted line ratios from each model 
are compared. The majo rity of the FIR fine structure line 
intensities are taken from lNegishi e t al. (200l]), who present 
ISO LWS observations for a large sample of galaxies. The 
line fluxes listed in table 2 of their paper have been converted 
to intensitie s using the extended source correction factors of 
iLlovdl (|2003l ) and an effective aperture of fib = 1.1 x 10 sr. 
For M51, the CO(1-0), CO( 3-2) and [Ci] 609 um integrated 
intensities listed in table 4 o f K ramcr et all (|2005h have been 
adopted, since these values are already at an appropriate 



angular resolution of 80 arcsec, derived by smoothing avail- 
able map data. The CO and [Ci] da ta for NGC6946 are 
taken from llsrael fc Baas! (|200ll . I2002T) and have been con- 
volved to the required res olution using an a ssumed source 
size of 26 arcsec x 23 arcsec (jNieten et al.lll999T ). The adopted 
line intensitie s for M82 were drawn f rom observ a tional data 
presented by IColbert et all (|l999l '). iMao et al.l (120001 ) and 
llsrael fc Baasl d2002l) for the FIR , CO and [C i] transitions, 
respectively. iBolatto et ail (|2000| ) have used line intensities 
observed toward nebula N27 of the Small Magellanic Cloud 
to estimate the physical conditio ns within the region by com- 
parison with the PDR models of lKaufman et all (|l999h . The 
ISO LWS data presented in their paper are used here, ad- 
justed to our adopted beam size of 75 arcsec, along with 
their observed value for the [C i] integrated intens ity in the 
region. The CO(1-0) map of iRubio et all l|l993l ) covers a 
smaller area within N 27, but is thought t o encompass most 
of the CO emission (|Bolatto et al.l I2000T ) and is therefore 
adopted here. Unfortunately, no CO (3-2 ) data are cu r rently 
available in the literature for this source. IRubio et alj (|l993r ) 
have mapped the CO (2-1) line in the region and we there- 
fore use these data to calculate the ratios [C ii]pdr/CO(2-1), 
[Ci] 609 u.m/CO(2-l) and CO(2-1)/CO(1-0) in place of the 
equivalent CO(3-2) line ratios. The CO and [Ci] data were 
convolved to the appropriate angular resolution assuming a 
source si ze of 110 arc s ec x 1 10 arcsec, determined from the 
maps of IRubio et~al] l|l993l '). The [Nil] 122 u.m data used 
to estimate the contribution of the i onized medium to th e 
[C ill 158 u.m e missio n ar e taken from [ Kramer et alj (|2005h . 
iNegishi et ail (|200ll ) and IColbert et all (|l999h forthe rele- 
vant regions in M51, NGC 6946 and M82, respectively. The 
resulting contributions from the ionized gas in these galaxies 
are found to range from ~30-40 per cent of the total [C n] 
emission, the greatest correction being applied to the M51 
data. No [Nil] data are available in the literature for the 
region of interest within the SMC. 

The line ratios predicted by the best fit model for each 
galaxy are listed in Tableland the physical parameters of 
each best fit model are listed in Table [S] Fig. [TJ{H show con- 
tours of reduced \ 2 , calculated for each grid model across the 
adopted range of densities and FUV fluxes, for each galaxy 
in our sample. White labelled contours mark the 1, 2 and 3<r 
confidence levels for each \ 2 distribution. The cloud age and 
size are fixed in each plot to the best fit values listed in Ta- 
ble [5] Models with Solar metallicity were initially used to fit 
the observed line ratios of M51, NGC 6946 and M82, which 
display similar metallicities. When determining the best fit 
models for the SMC region, the analysis was restricted to 
only those models with Z/Zq = 0.1, appropriate for this 
low metallicity environment. 



4.1 M51 

For the central region of M51, in which our adopted beam 
size encompasses an area approximately 3.5 kpc across, 
we find that the observed line intensity ratios are best fit 
by average conditions characterised by reasonably dense 
(n = 7 x 10 3 cm -3 ) gas immersed in an ambient FUV field 
(Go = 20 Habing). T hese values agr e e rem arkably well with 
those determined by iKramer et all l|2005h . who perform a 
similar \ 2 fit to th e observed line ratios using the PDR 
models presented bv lKaufman et all (|l999h . They find that 
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Table 3. Observed line ratios for each galaxy in the sample. 



Galaxy 


LO ij 63 
[Cn| pnn 


[Oi] 145 
[OiJ 63 


[CiJ 609 


I^WJPDR 

CO(3-2) 


[Cl] 609 
CO (3-2) 


CO(3-2) 
CO(1-0) 


M51 


0.730 


n/a 


39.4 


43.9 


1.114 


14.74 


NGC 6946 


0.935 


0.058 


112.2 


58.2 


0.519 


14.11 


M82 


1.329 


0.092 


128.2 


605.6 


4.726 


16.11 


SMC N27 


0.873 


0.026 


219.7 


512.3* 


2.331* 


8.694* 



*CO(3-2) data were not available at the location of the N27 nebula so these line 
ratios were instead c alculated using the CO(2— 1) integrated intensity derived from the 
observational data of lRubio et alj l|l993h . 



Table 4. Line ratios predicted by the best fit model for each galaxy in the sample. 

T~~~ 1 1 1 63 Oil 145 [Cii| PDR [Cii| PDR |Ci| 609 CO(3-2) 

_ 7 [Cirj PTm [OiJ 63 ]Ci] 609 CO(3-2) CO(3-2) CO(l-O) 

M51 0.728 0.037 49.7 54.0 1.087 21.84 

NGC 6946 0.885 0.036 49.2 45.7 0.928 23.42 

M82 2.081 0.042 85.4 500.1 5.854 14.03 

SMC N27 0.488 0.021 209.4 1113* 5.315* 6.218* 

*Calculated using the CO(2— 1) integrated line intensity to allow comparison with the 
observed line ratios. 
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Figure 1. Centre of M51. Contours of the reduced \ 2 fit to the 
data are shown across the range of parameter space covered by 
the grid of models in the density-radiation (n— Go) plane. The 
cloud age and size are fixed at the best fit values for M51 listed 
in Table IBI Contour levels are 2.2 (black), 2.7, 4.0 (lrr confidence 
level), 7.8 (2a), 13.5 (3<r), 25, 50, 100, 200, 500, 1000 and >1000 
(white) . 



local densities of n ~ 10 4 cm~ 3 and FUV fluxes of Go ~ 20 
Habing best fit the observations, both at the central position 
and within the spiral arms of the galaxy. This level of agree- 
ment is encouraging, suggesting that the method is reliable 
and that our models are consistent with those of others. 
However, our approach has the additional benefit of deter- 
mining the average cloud size and estimated degree of chem- 
ical evolution within the region, since the models produce 
line intensities that are both depth- and time-dependent. We 
find that the line ratios are best fit by cloud models with vi- 
sual extinctions of Ay ~ 10 mag and ages of ~10 Myr (when 
evolving from initially atomic gas). Such conditions are ex- 
pected in dense regions where many clouds along the line 
of sight give rise to optically thick emission (with the pos- 
sible exception of the [Cl] and [Cn] lines) and the gas is 



maintained in molecular form. Our best fit model predicts 
a cosmic-ray ionization rate of £ = 5xl0~ 17 s _1 , consistent 
with the low radiation field and q uiescent star fo r matio n ac- 
tivity in this region. The fact that lKramer et al | i|2005l ) find 
similar conditions in the spiral arms of M51 suggests that 
our models might be equally applicable to these regions. 

The centre of M5 1 is known to have a metallicity slightly 
above Solar (see, e.g.. lZaritskv et al]|l994l ). To see if a better 
fit to the observed line ratios might be obtained by using a 
higher metallicity, we computed additional models with the 
same density, incident FUV flux and cosmic-ray ionization 
rate as the best fit model, but with metallicities of Z = 2 Z@, 
3 Zq , 4 Zq and 5 Zq (the original best fit model was com- 
puted at Solar metallicity). The resulting reduced \ 2 val- 
ues for these models are 1.62, 1.63, 1.73 and 1.74 for the 
Z — 2 Zq to 5 Zq cases, respectively. Models with Z = 2 Zq 
and Z = 3 Zq show a slight improvement in the goodness 
of fit, with reduced x 2 values lower than that of the original 
best fit model (Xmin — 1.68), but given the uncertainties in 
the line strengths, these differences are unlikely to be signifi- 
cant. It is probable that these line ratios are only able to pro- 
vide loose constraints on the metallicity. Since the Z = 2 Zq 
model produces the best fit, we adopt this metallicity for the 
central region of M51, however, we note that, in the absence 
of additional constraints, the uncertainty in the metallicity 
will lead to an uncertainty in the derived value of Xco- The 
minimum reduced \ 2 for this model i s a fac tor of 2 better 
than that achieved by iKramer et ail (|2005T ). We therefore 
regard this model as representative of the conditions within 
interstellar molecular gas in normal spiral galaxies. 



4.2 NGC 6946 

Similar physical parameters are found to best fit the line 
ratios determined for the central region of NGC 6946 (n — 
10 4 cm" 3 , Go = 20 Habing, t ~ 10 Myr, A v ~ 10 mag 
and £ = 5xl0 -17 s _1 ). This is surprising, given that the 
galaxy possesses a starburst nucleus and higher FUV fluxes 
would therefore be expected. We note that the fit is not 
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Table 5. Best fit physical parameters for each of the four galaxies, based on the model with 
the lowest \ 2 value. 



Galaxy 


n (cm 3 ) 


Go (Habing) 


* (yr) 


A v (mag) 


Z/Z Q 




X 2 - 

v mm 


M51 


7+3 1 X 10 3 


20+J2 


10 v 


9.8 


2 


5xlO- iY 


1.6 


NGC 6946 


lto.3 >< 104 


20±H 


10 7 


9.9 


3 


5xicr 17 


6.0 


M82 


lt° - 9 i x 10 3 


500±ii 


10 s 


5.3 


1 


5xl0~ 17 


2.9 


SMC N27 


Ho.4 x 10 3 


lOOljjf 


10 s 


2.8 


0.1 


5xl0~ 17 


1.9 
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Figure 2. Centre of NGC 6946. Contours of the reduced \ 2 fit to 
the data are shown across the range of parameter space covered 
by the grid of models in the density-radiation (n-Go) plane. The 
cloud age and size are fixed at the best fit values for NGC 6946 
listed in Table[S] Contour levels are 7.9 (black), 8.4, 9.7 (la confi- 
dence level), 13.6 (2a), 19.2 (3<x), 25, 50, 100, 200, 500, 1000 and 
> 1000 (white). 

as good here as it was for M51, with a minimum reduced 
X of 7.41 for the best fit model and some degree of uncer- 
tainty in the appropriate value of Go, as demonstrated by 
the elongated contours in Fig. [2] Since the m etallicity in the 
centr al region of NGC 6946 is above Solar jZaritskv et al.l 
1 19941 ). we consider models with higher metallicities in the 
same manner as for M51. We find that increasing the metal- 
licity to Z = 3 Zq produces a better fit (xmin = 5.99), whilst 
higher metallicities yield poorer fits (\ 2 > 6). We therefore 
adopt Z = 3 Zq as the best fit model for the central re- 
gion of NGC 6946. Again, we note that the metallicity is 
only weakly constrained by this technique and the choice of 
app ropriate value is somew hat uncertain. 

IContursi etaD (j2002h fitted the line intensity ratio 
[C ii]pdr/[0 i] 63 [im, together with the line-to-FIR contin- 
uum ratio, ([C ii]pdr+[0 i] 63 urn) /FIR, for ISO IWS obser- 
vations of various positions within NGC 6946, aga i n emp loy- 
ing the PDR model predictions of lKaufman et all (|l999l ). At 
the central position, they derived for the neutral gas an av- 
erage density of n fs 2 x 10 3 cm -3 and an incident FUV 
flux of Go ~ 900 Habing. These values are significantly dif- 
ferent to those that we have derived, predicting the gas to 
be a factor of 5 lower in density and, more dramatically, the 
incident FUV flux to be close to a factor of 50 higher. Such 
a large difference suggests that, by comparing different line 
ratios, we may be probing different gas components within 
the region. The large beam size, corresponding to an area 
~2 kpc across at the distance of NGC 6946, makes it likely 
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n [cm 3 ] 

Figure 3. Centre of M82. Contours of the reduced \ 2 fit to the 
data are shown across the range of parameter space covered by 
the grid of models in the density-radiation (n— Go) plane. The 
cloud age and size are fixed at the best fit values for M82 listed 
in Table[3] Contour levels are 3.4 (black), 3.9, 5.2 (la confidence 
level), 9.1 (2a), 14.7 (3a), 20, 25, 30, 35, 42, 50 and > 50 (white). 

that the emission arises from a mixture of active and less 
active regions. The FIR continuum traces the strength of 
the FUV flux, which is absorbed by the dust and re-emitted 
in the IR. However, whilst this emission is optically thin, 
and therefore seen across the whole region, the FUV flux 
responsible is rapidly attenuated in the clouds it is incident 
upon, so the surface area of the gas that is actually sub- 
jected to the strong FUV flux may constitute only a small 
fraction of the total area. By using this diagnostic to derive 
the physical conditions, the inferred cloud properties may 
therefore be biased towards those of only a small fraction 
of the total gas content. Both sets of physical parameters 
are likely to represent the conditions within certain regions 
of the galaxy and we therefore consider both when deriving 
conversion factors in the next section. 

4.3 M82 

The starburst galaxy M82 is found to have average gas prop- 
erties of n — 10 3 cm -3 and Go = 500 Habing within its cen- 
tral region (an area ~1 kpc across is enclosed in the 75 arcsec 
beam at this distance), according to the best fit model de- 
termined by our \ 2 analysis (Xmin = 2.9). These pr o pertie s 
are consistent with those derived by IColbert et al.l l|l999h , 
who find n » 2000 cm~ 3 and Go ~ 600 Habing using a 
combination of H II and PDR models to self-consistently fit 
the observed FIR ionic and atomic fine structure lines. A 
typical cloud age of ~100 Myr is predicted by the best fit 
model. We find that the observed line ratios are best fit by 
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Figure 4. N27 region of the Small Magellanic Cloud. Contours 
of the reduced \ 2 fit to the data are shown across the range of 
parameter space covered by the grid of models in the density- 
radiation (n-Go) plane. The cloud age and size are fixed at the 
best fit values for the region listed in Table [5] Contour levels are 
5.4 (black), 5.9, 7.2 (ler confidence level), 11.1 (2a), 16.7 (3a), 20, 
30, 40, 60, 80, 100 and > 100 (white). 

smaller clouds, with total visual extinctions of Av ~ 5 mag. 
This is in ag r eemen t with the recent observational study of 
iFuente et alj |2005), who find that the observed abundance 
ratio [CN]/[HCN] > 5, determined for the central 650 pc 
of M82, requires cloud sizes of Av < 5-6 mag. For com- 
parison, we plot the [CN]/[HCN] abundance ratio produced 
by our best fit model as a function of visual extinction in 
Fig. [5] From this plot, it can be seen that a cloud size of 
Av ;$ 4 mag is also required to maintain an abundance ra- 
tio > 5 in this model. 

ISuchkov, Allen fc Heckmanl |l993) proposed an in- 
creased cosmic-ray ionization rate of ( = 4 x 10 -15 s _1 
in the centre of M82 to account for the observed physical 
conditions of the molecular gas. The grid models with the 
same density and incident FUV flux as the best fit model, 
but with cosmic-ray ionization rates of f = 5 X 10~ 16 and 
5 x 10~ 15 s _1 (i.e. 10 and 100 times the best fit rate) pro- 
duce reduced x 2 values of 8.7 and 5.5, respectively. Both are 
higher than the best fit model value of Xmm = 2.9, although 
the model with a rate 100 times higher shows a better fit 
than the model with only a factor 10 increase in the cosmic- 
ray ionization rate. Since no improvement in the goodness 
of fit is obtained by using a higher value of £, we adopt 
the physical parameters of the original best fit model when 
determining appropriate conversion factors for starburst re- 
gions. 



10 2 




10 1 

2 4 6 8 10 

A v [mag] 

Figure 5. Centre of M82. The abundance ratio of CN to HCN 
([CN]/[HCN]; solid line) as a function of visual extinction Av for 
the best fit model (n = 10 3 cm" 3 , Gp = 50 Habing), compared 
to the ratio observed bv lFuente et al. | j2005l , dotted line). 

sonably dense gas (n = 10 3 cm~ 3 ) subjected to elevated 
FUV fluxes (Go = 100 Habing), likely due to the active 
star formation and reduced extinction in this low metallic- 
ity environment. Thes e results are in very good agreement 
with those obtained bv lBolatto et al.l (|2000l ) . determined us- 
ing FIR fine structure lines and continuum emission to be 
n ~ 300-1000 cm -3 and Go ~ 30-100 Habing. Furthermore, 
whilst previous models have had to impose the limit of low 
Av observed in these clouds, this property is independently 
predicted as a result of our fitting procedure (Av ~ 2.8 mag 
in our best fit model). The line ratios observed in N27 are 
less well fit than those of M51 and M82 (xLn = 4.9 for the 
best fit model of N27) . An analysis of the agreement between 
model predictions and observations shows that the predicted 
[Cn]/CO(2-l) ratio is the most discrepant of the six line ra- 
tios considered, being over a factor of 2 higher than the ob- 
served value. We note that in our adopted line ratios for N27 
we have not accounted for the ionized gas contribution to the 
[Cll] 158 nm emission, since the required [Nil] 122 u.m data 
are unavailable in the literature. However, subtracting the 
ionized gas contribution would decrease the observed value 
of the [C ll]/CO(2-l) ratio, which would actually reduce the 
quality of the fit. Instead, a more favourable explanation for 
the disagreement between predicted and observed line ratios 
is the high O/C abundance ratio in the SMC, a factor 2.5 
greater than the local value (jWilson et al.| [2005). This was 
not taken into account in the models. We adopt the physi- 
cal parameters of the best fit model for our consideration of 
suitable conversion factors for dwarf irregular galaxies. 



4.4 SMC (N27) 

In the low metallicity environment of the Small Magellanic 
Cloud, the molecular gas is known to exhibit significantly 
different properties to those typically found in regions with 
metallicities closer to Solar. By restricting our parameter 
search to models with metallicity Z = 0.1 Zq (the aver- 
age value derived for the SMC by iDufourl Il984l ) we obtain 
a best fit model for the line ratios observed in the SMC 
nebula N27 that predicts the presence of clouds with rea- 



4.5 Discussion of best fit parameters and their 
uncertainties 

The best fit ages determined for each galaxy provide an in- 
dication of the degree of chemical evolution from atomic gas. 
Whilst ages of 10 7 -10 8 yr are predicted for M51, NGC6946 
and M82, suggesting that the emission arises from predom- 
inantly molecular gas, a relatively young age of 10 5 yr is 
derived for the SMC, which implies that the gas in this re- 
gion is mainly atomic, either due to the slower formation 
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rate of H2 in the low metallicity environment or because 
the reduced shielding results in greater photodissociation 
of molecular species. The radiation field is stronger still in 
M82, yet the molecular gas is predicted to be well-evolved. 
This is probably because M82 has a higher metallicity and 
a higher typical Ay than are found in the SMC, and both 
of these promote chemical evolution. The degree of chemical 
evolution can have a strong effect on the CO-to-Hb conver- 
sion fac tor, reducing the value of Xco at early times (see 
I Paper J . for details). The cosmic-ray ionization rate for each 
best fit model is found to take its standard value in all cases, 
suggesting that the line ratios considered here may not be 
sensitive to this parameter. Whilst the \ 2 analysis indicates 
that models with metallicities above Solar may be more ap- 
propriate for the centres of M51 and NGC 6946, these metal- 
licities cannot be tightly constrained with our choice of line 
ratios, as discussed in previous sections. 

To provide some estimate of the uncertainty in the best 
fit density and radiation field strength for each galaxy, 1, 2 
and 3<r confidence level contours are overlaid on the plots 
of reduced % 2 m Fig- HHS an d the range of each parameter 
contained in the la region is taken as the uncertainty in 
that parameter, as listed in Table [5] The uncertainty in the 
best fit gas density is typically found to be a factor of 2- 
3, although for region N27 in the SMC the density is less 
well constrained and may be up to a factor of 5 higher. 
The estimated FUX fluxes suffer a slightly greater degree of 
uncertainty, generally being constrained to within a factor 
of 4. Again, the greatest degree of uncertainty in Go is found 
for the SMC, consistent with the relatively poor fit that was 
obtained for this galaxy. Given the weak constraints on the 
metallicity provided by our fitting procedure, we also adopt 
an uncertainty of a factor of 2 in metallicity. The estimated 
uncertainties in these parameters will be used in Section [6] 
to examine the variation in derived conversion factors. 

In the next section we use the best fit models for each 
galaxy in our sample to derive CO-to-H2 and C-to-Hb con- 
version factors that are appropriate for molecular gas within 
these environments and suggest transitions of CO that might 
yield the most reliable mass estimate in each case. 



5 APPROPRIATE CONVERSION FACTORS 
FOR MOLECULAR MASS ESTIMATES IN 
DIFFERENT GALAXY TYPES 

Having determined the physical conditions that best repre- 
sent the average properties of the neutral ISM within each 
galaxy type, we now derive CO-to-H2 and C-to-Hb conver- 
sion factors appropri ate for those conditions in the manner 
described in iPaper it In addition to the standard value of 
Xco, relevant to CO J = I — > emission, we extend our 
study to include conversion factors for higher rotational lines 
of CO and for the [C i] 609 [im line that might serve as use- 
ful alternatives under certain conditions. The UCL_pdr code 
calculates the integrated intensities of all CO rotational lines 
up to J = If — ^ 10 and values of Xco can therefore be de- 
rived for each transition using a more general definition of 
the CO-to-Eb ratio, 
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Figure 6. Normal spiral galaxies. Xco versus Ay profiles for 
the best fit physical parameters for the centre of M51. The val- 
ues of Xqo shown are for emission from the various rotational 
transitions of CO. 



where jV(H2) is the column density of H2, as before, and 
Ta(CO J = u — > V) is the antenna temperature of the CO 
J = u — » I line. We will use the notation Xqq 1 to distinguish 
between the different CO-to-H2 conversion factors through- 
out this section. The C-to-H2 conversion factor Xc 1 for each 
galaxy can be derived in a similar way, using the [C 1] 609 u.m 
integrated intensities calculated by the ucl_pdr code. 

Fig. IBl HOI show Xco as a function of visual extinction 
Ay for various transitions of CO, from the best fit model for 
each galaxy in our sample. Recall that, for these Xco pro- 
files, the observationally-relevant conversion factor in each 
transition corresponds to the profile minimum, which rep- 
resents the peak em ission from the cloud (or ensemble of 
clouds; see IPaper ill . Fig. [TT] shows Xc 1 as a function of 
Ay for each galaxy. Again, the profile minima represent 
the conversion factors relevant to observations. Table [5] lists 
the minimum values of Xco and Xci [in units of 10 20 
cm~ 2 (Kkms -1 ) -1 ] for each galaxy type in various tran- 
sitions of CO (1-0, 2-1, 3-2, 4-3, 6-5 and 9-8), together 
with the cloud depths at which they occur (in Ay). These 
are therefore the appropriate values of Xco and Xci that 
should be used to derive molecular mass estimates in these 
regions. The results for each of the four galaxy types are 
discussed briefly below. For the remainder of this section, 
the units of the CO-to-H2 and C-to-Eb conversion factors, 
cm" 2 (Kkms -1 )" 1 , are omitted to save space. 



5.1 Normal spiral galaxies 

We adopted the conditions within the central ~3 kpc of 
M51 as characteristic of those within normal late-type spi- 
ral galaxies, both in their centres and in their spiral arms. 
Using the best fit model for these conditions, we find a value 
for the CO(1-0) conversion factor of X C q = 2.5 x 10 19 . 
This value is an order of magnitude lower than the canoni- 
cal CO-to-Eh conversion facto r derived for the S olar neigh- 
bourhood, X C q° = 2 x 10 20 (|Dame et al.ll200lh . The drop 
in Xco values is due to both the increased density and 
high er meta llicity, which result in stronger CO emission 
(see IPaper J . for details). Observational studies of the cen- 
tres of nearby spiral galaxies have consistently inferred val- 
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Figure 7. Weak starbust nuclei. Xqq versus Ay profiles for the 
best fit physical parameters for the centre of NGC6946. Values 
of Xqq for the various rotational transitions of CO are shown. 




- 7 = 4^3 : 



: 7-2 >1 ./-6 >5 : 

ln i 9 7=3^2 7=9^8 " 

1U I i i i i i i i i i i i i i i i i i i i 

2 4 6 8 10 

A v [mag] 

Figure 8. Starburst nucle i. Xqq versus Ay pro files for the phys- 
ical parameters derived bv lContursi et al. (2002) for the centre of 
NGC6946. Values of Xqq for the various rotational transitions 
of CO are shown. 



ues of Xm° that are lower than the canonical value (see, 
e.g..|Malonev fc Blacklll98Sl : ISorai et alJkood : llsrael fc Baasl 
12003 , and references there in), in agreement w ith our theo- 
retical value. Furthermore, llsrael et alj l|2006l ) have recently 
determined the CO-to-Hb conversion factor for the centre of 
M51 to be Xqq = 5.0 ± 2.5 x 10 19 from observations of 
neutral carbon and CO emission within the region. This is 
consistent with the value obtained from our best fit model 
and we feel that this independent observational verification 
shows our theoretical approach to be reliable. Looking to 
the higher transition lines of CO, the CO (2-1), CO (3-2) and 
CO(4-3) lines all give similar conversion factors (2.0xl0 19 , 
2.4xl0 19 and 3.7xl0 19 , respectively) and become optically 
thick (r > 1) beyond A v « 2 mag. The CO (6-5) and CO (9- 
8) lines remain optically thin to greater cloud depths, but 
are fairly weak due to the low gas temperatures (~150-10 K 
from the surface to the cloud interior) . The lower transitions 
are therefore more likely to yield reliable molecular mass es- 
timates, particularly when several lines intensities are used 
in combination. The appropriate C-to-H2 conversion factor, 
determined from the best fit model, is Xq, = 5.9 x 10 19 
and its variation with visual extinction displays a similar 
trend to that of the CO(1-0) conversion factor, though the 
minimum occurs nearer the cloud surface (see Fig. [XT]) . 



5.2 Starburst nuclei 

Our best fit physical parameters for the centre of NGC 6946 
are similar to those obtained for M51, representative of more 
quiescent interstellar gas, although a slightly higher density 
(10 4 cm -3 instead of 7 x 10 3 cm" 3 ) and higher metallic- 
ity (3 Zq instead of 2 Zq) axe predicted. As such, the de- 
rived CO-to-Hb conversion factors are also similar, though 
they are typically 20-50 per cent lower than those found 
for M51 (see Fig. [7] and Table©. The best fit model pro- 
duces a CO (1-0) conversion factor X C q° = 1. 6 x 10 19 and a 
C-to-H 2 conversion factor Xci = 3.5 x 10 19 . llsrael fc Baas! 
|200ll ) derived a value of X C q° ~ 1 x 10 19 for the centre 
of NGC 6946, which is remarkably similar to our theoretical 
value. The depth-dependent profile of Xq i is similar to that 



of Xco for the lower CO transitions, as was the case for 
M51. 

As discussed in the previous section, these conditions 
may be selectively sampling a particular component within 
the galaxy centre and may not be fully representative of typ- 
ical nuclear starburst activity. We therefore also consider an 
alternative model for the starburst nucleus of NGC 6946, 
based on the parameters determined by IContursi et al.l 
||2002D , i.e. n = 2 x 10 3 cm" 3 and Go = 10 3 Habing. The 
resulting Xco profiles for the various CO transitions are 
shown in Fig. [8] and the appropriate conversion factors to 
be used with each transition to derive molecular masses are 
listed in Table [6] The increased FUV radia tion striking the 
interstellar gas in the lContursi et al.1 (|2002h model leads to 
enhanced photodissociation of CO, pushing the C + /C/CO 
transition deeper into the cloud, where the lower gas tem- 
peratures (~40 K) result in reduced emission strengths. The 
CO(1-0) conversion factor is Xqq = 2.0 x 10 20 , the value 
obtained for local regions within the Milky Way. From exam- 
ining Fig.[8j it can be seen that the CO (4-3) line may prove 
to be the best tracer of total molecular mass in this environ- 
ment, since the Xqq 3 conversion factor remains fairly con- 
stant to high visual extinction. The appropriate value for the 
conversion factor in this case is Xqq 3 =1.0 x 10 21 . The cor- 
respon ding C-to-Eb conversion factor for the lContursi et aD 
(2002) model is X Cl = 1.2 x 10 20 . Its depth-dependent pro- 
file (shown in Fig. Illfl displays a more pronounced minimum 
than that of the CO (4-3) line, making it slightly less reliable 
as a tracer of total molecular mass under these conditions. 



5.3 Irregular starburst galaxies 

ISmith et al.l (|l99ll ) have derived the CO-to-H2 conversion 
factor for the central region of M82 from observations of the 
dust continuum. They find a value of ~ 1.2 x 10 20 , which is 
in good agreement with the value obtained from our best 
fit model, X C q° = 1.5 x 10 20 . The X co profile minima for 
the low transitions of CO occur deep within the cloud at 
Av ~ 6 mag (see Fig. |9j, as in the previous case, due to 
the higher FUV flux (500 Habing). The conversion factors 
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- 7 = 4^3 : 



: 7-2 >1 ./-6 >5 : 
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1U I i i i i i i i i i i i i i i i i i i i 

2 4 6 8 10 

A v [mag] 

Figure 9. Irregular starburst galaxies. Xqo versus Ay profiles 
for the best fit physical parameters for the centre of M82. Values 
of Xco f° r the various rotational transitions of CO are shown. 




- 7 = 4^3 : 
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2 4 6 8 10 
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Figure 10. Dwarf irregular galaxies. Xqo versus Ay profiles 
for the best fit physical parameters for the N27 nebula of the 
Small Magellanic Cloud. Values of Xqq for the various rotational 
transitions of CO are shown. 



for transition lines up to CO (4-3) all show fairly constant 
values for visual extinctions above ~5 mag, suggesting that 
accurate molecular masses may be obtained for larger cloud 
sizes. However, the best fit model predicts typical cloud sizes 
of Ay < 5 mag for M82, in accordance with observations, 
so the appropriate CO-to-H2 conversion factors for clouds 
of low Ay are likely to vary quite significantly, as shown by 
the sharp rise in Xco towards the left hand side of Fig. [9] 
The difficulties of using the CO-to-Hb conversion factor in 
suc h translucent clouds h ave been discussed extensively (see, 
e.g.. lMagnani et alj|2003l . and references therein). The C-to- 
H2 conversion factor, however, is less variable at low visual 
extinction (see Fig. II 1|) and its minimum value, Xci = 1.3 x 
10 20 , occurs at Ay = 4.3 mag, making it a better tracer of 
the molecular mass in these clouds. 



5.4 Dwarf irregular galaxies 

The best fit model for the N27 nebula within the Small 
Magellanic Cloud can provide conversion factors suitable for 
metal-poor dwarf irregular galaxies. It predicts a CO(l-O) 
conversion factor of X C o° = 9.5 x 10 20 , approximately 5 
times greater than the canonical value for the Milky Way. 
This is primari ly due to the low metallicity (Z = 0.1 Zq), 
as discussed in iPaper j . The low visual extinction in these 
clouds may mean that the appropriate value is higher still, 
since Xco increas es nearer the cloud surface (see Fig. I10[) . 
iLerov et all (|2006l ) derive a CO-to-H 2 conversion factor of 
X C q° = 6± 1 x 10 21 for regions within the SMC that are oc- 
cupied by CO, using new FIR continuum maps to determine 
the dust mass. This is a factor of 6 higher than the theoreti- 
cal value predicted by our best fit model and may imply that 
the treatment of metallicity-dependence in the model does 
not fully account for the reduction in CO emission needed to 
raise the value of Xco (metallicity is the physical parameter 
that cau ses the greatest variation in the value of Xco; see 
I Paper J) . Alternatively, if typical cloud sizes in the SMC are 
Ay ~ 1.6 mag, rather than the 2.8 mag predicted by the 
best fit model, then the appropriate CO-to-Hb conversion 
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SMC N27 




10' 



10 



Ay [mag] 



factor would be X C o ° = 6 x 10 21 , due to the steep increase 



Figure 11. Xrji versus Ay profiles for each galaxy. The profile 
labelled NGC 6946a is produced by the best fit physical parame- 
ters for the centre of NGC 6946 obtained in this work; the profile 
labelled NGC 69466 is produc ed by the physical parameters de- 
rived by IContursi et al.1 d2002h . 



in Xco at l ow A y. Such cloud sizes are indeed found by 
ILerov et al.1 (|2006l ). 

The C-to-H2 conversion factor derived from the best fit 
model is Xci = 1-9 x 10 20 , however, the minimum of the 
Xc 1 profile occurs at Ay ~ 4 mag (see Fig. II 1 p , deeper into 
the cloud than that of the CO-to-Hb profiles, making it less 
appropriate as a tracer of molecular mass in this case. 

The conversion factors listed in Table [6] provide a useful 
reference for observers, presenting 'standard' conversion fac- 
tors that are more appropriate to these diverse environments 
than the canonical value derived for the Milky Way. By ob- 
taining data in multiple transitions of CO and C I, several 
conversion factors can be used to provide better constraints 
on the estimated molecular mass within a region. 
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Table 6. Recommended values for the CO-to-H2 and C-to-H2 conversion factors for each galaxy type. CO-to-H2 conversion factors are 



listed for various transitions of CO. 


Values of Xco and Xq 


l are given 


in units of 10 20 


cm 2 


(Kkms" 1 )- 


■ l 








Galaxy type 


J = 


1 -s-0 


J = 


2 -> 1 


J = 


3^2 


J = 


4^3 


J = 


6 5 


J = 


9^8 


[Ci] 


609 \aa 




A v 


Xca 


A v 


^CO 


A v 


^CO 


A v 


■X"co 


A v 


Xco 


A v 


Xco 


A v 


Xci 


Normal spiral 


2.6 


0.3 


2.3 


0.2 


2.2 


0.2 


2.3 


0.4 


2.4 


1.2 


1.8 


4.2 


1.8 


0.6 


Starburst nucleus" 


2.4 


0.2 


2.2 


0.1 


2.1 


0.1 


2.1 


0.2 


2.3 


0.6 


1.6 


2.3 


1.7 


0.4 


Starburst nucleus 6 


7.4 


2.0 


6.7 


2.1 


5.5 


3.5 


5.2 


10.4 


3.6 


25.9 


3.5 


61.2 


4.8 


1.2 


Irregular starburst 


5.9 


1.5 


5.5 


1.7 


5.5 


2.9 


5.7 


9.5 


4.1 


79.4 


4.1 


203.7 


4.3 


1.3 


Dwarf irregular 


2.1 


9.5 


2.1 


12.0 


2.3 


37.7 


2.3 


179.3 


2.2 


203.9 


2.2 


508.0 


3.S 


1.9 



a Determined for the best fit model obtained for NGC 6 946 in this work (see S ection [4]l. 
b Model values calculated for the physical parameters of lContursi et ah (2002). 



6 DISCUSSION AND CONCLUSIONS 

We have determined appropriate physical parameters that 
represent the average conditions within four key extragalac- 
tic environments; the interstellar gas within the spiral arms 
and nuclei of normal late-type spiral galaxies (modelled on 
the central region of M51), the starburst nuclei of spiral 
galaxies and more extreme starburst environments (includ- 
ing the centre of the barred spiral galaxy NGC 6946 and the 
prototypical starburst galaxy M82) and the low metallici- 
ties and high UV fields typical of dwarf irregular galaxies 
(in this case, modelled on region N27 of the Small Magel- 
lanic Cloud). Best fit models for these environments have 
been determined by fitting observed line intensity ratios in 
the selected regions for a range of FIR and sub-mm diag- 
nostic lines. 

The resulting models have been used to derive appro- 
priate CO-to-Eb and C-to-H~2 conversion factors that will 
allow more reliable molecular mass estimates to be inferred 
from observations of CO rotational lines and the [C i] 609 \im 
fine structure line in a wide variety of extragalactic environ- 
ments. These conversion factors are listed in Table H3 Con- 
version factors are presented for a range of CO transitions, 
enabling estimates of the molecular mass to be further con- 
strained using several line intensities simultaneously. Fur- 
thermore, different emission lines may probe specific com- 
ponents within the interstellar medium, e.g. dense and hot 
gas in the case of the high J CO transitions. 

Uncertainties in the best fit physical parameters have 
been determined for each galaxy, and these are now used to 
examine the corresponding variation in the derived CO-to- 
H2 conversion factors. The uncertainties in the best fit gas 
density and FUV flux are listed in Table [5] and we assume 
a factor of 2 uncertainty in the metallicity. Uncertainties 
in the other physical parameters are not considered here. 
Xco decreases with increasing density, as the CO emission 
becomes stronger nearer the cloud edge, whilst it generally 
increases with FUV flux, as CO is more readily photodis- 
sociated; Xco depends strongly on metallicity, decreasing 
as metallicity increases and CO emission be comes str onger 
(for a detailed discussion of these trends, see lPaper J ). The 
sensitivity of the C-to-H"2 conversion factors to these uncer- 
tainties is not discussed here, but is found to be very similar. 

For the centre of M51, the uncertainty in the gas density 
typically causes Xco to vary by less than a factor of 2, with 
the exception of the higher CO transitions (J = 6 — > 5 and 
above), for which variation by over a factor of 3 is found. 
Uncertainty in the best fit FUV flux has little or no influ- 
ence on Xco, with at most a factor of 2 variation in the 



J = 6 — » 5 and 9^8 conversion factors. The assumed 
uncertainty of a factor of 2 in metallicity has the greatest 
impact on the derived values of Xco , causing them to vary 
by up to a factor of 4. The higher transition lines of CO are 
more sensitive to variations in the gas density and temper- 
ature, since their critical densities for excitation are higher 
and greater temperatures are required to excite them. For 
the best fit parameters for NGC 6946, the uncertainty in 
the density leads to variations in Xco of up to a factor of 
2, though without the increase in the higher transition lines 
that was found for M51. The slightly larger FUV flux uncer- 
tainty leads to a marginal increase in the degree of variation 
of Xco, though it is still typically no more than a factor 
of 2. Variation in the metallicity causes up to a factor of 4 
variation in Xco- 

Relatively tight constraints on the best fit density and 
radiation field strength for M82 mean that the resulting vari- 
ation in Xco is also small and always less than a factor of 2. 
Assuming that the metallicity is again uncertain by a factor 
of 2, larger variations in Xco are produced than were seen 
in M51 or NGC 6946 (X C o varies by a factor of 3-4 in the 
lower transitions and by over 6 in the 6 — * 5 and 9 — * 8 tran- 
sitions). The larger variation with metallicity seen in M82 
is probably due to the elevated FUV flux, which, combined 
with the change in shielding that accompanies the change 
in metallicity, can enhance the photodissociation of CO. 

The best fit physical parameters determined for the 
SMC suffer the greatest uncertainties and this is reflected 
in the degree of variation in the conversion factors. The low 
metallicity adopted for this galaxy (Z — 0.1 Zq) also causes 
Xco to be more sensitive to change in the physical param- 
eters. Therefore, the large uncertainty in the density (up to 
a factor of 5) leads to variations in Xco that range from a 
factor of 4 for the J = 1 — » transition to over a factor of 
35 for the 9 — > 8 transition line. The weak dependence of 
the CO-to-H2 conversion factor on the FUV flux means that 
the uncertainty in this parameter produces only a factor of 2 
variation in Xco ■ Uncertainty in the metallicity causes Xco 
to vary by a factor of 4-5. 

Overall, the uncertainty in the best fit density and radi- 
ation field strength for each galaxy typically cause the corre- 
sponding CO-to-Eb and C-to-Eb conversion factors to vary 
by no more than a factor of 2-3, whilst the adopted uncer- 
tainty in the metallicity has a stronger influence, producing 
a factor of 4 or higher variation in Xco and Xci- Conver- 
sion factors for the higher transition lines of CO are found to 
be more susceptible to the uncertainties in the physical pa- 
rameters. The conversion factors derived for the SMC suffer 



© 2007 RAS, MNRAS OOO.ITlfiH 



Conversion factors for specific galaxy types 13 



the greatest level of uncertainty, a consequence of the poorer 
quality of fit that was obtained for this galaxy. 

The essence of the X-factor method is that a single mea- 
surement of CO intensity - representing a wide range of 
physical conditions in a region of space - can be used to 
estimate the mass of gas in that region. Our method pre- 
dicts the most probable values for the physical conditions 
within the region considered, and hence the most probable 
value of the CO (or other) emission. We believe, therefore, 
that our method is consistent with the X-factor method and 
that our derived conversion factor values are appropriate to 
use, rather than other values from the possible range of Aco 
obtained as above. 

In future work, the range of physical parameters over 
which each specified conversion factor remains valid will be 
examined in more detail, providing limits on their applica- 
bility. The emission lines of other species might also prove to 
be better tracers of the molecular mass under certain con- 
ditions (e.g. HCN for dense gas); this prospect will also be 
considered in a future study. It is hoped that, taken together, 
these conversion factors will provide a firm basis upon which 
increasingly reliable molecular mass estimates can be made. 
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